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EIGENVALUE ASYMPTOTICS, INVERSE PROBLEMS 

AND A TRACE FORMULA FOR THE LINEAR 

q^-, DAMPED WAVE EQUATION 

O 
O 

<N 

>>: 

C$ • Abstract. We determine the general form of the asymptotics for 

Dirichlet eigenvalues of the one-dimensional linear damped wave 
operator. As a consequence, we obtain that given a spectrum corre- 
sponding to a constant damping term this determines the damping 
term in a unique fashion. We also derive a trace formula for this 
problem. 



Oh 

1. Introduction 

Consider the one-dimensional linear damped wave equation on the 
interval (0, 1), that is, 

> \ 

^ . ( w tt + 2a(x)w t = w xx + b(x)w, x e (0, 1), t > 

<N : (l.i) I w(0,t) = w(l,t) = 0, t>0 

■ 

^ I l w(x, 0) = w (x), w t {x,0) =w 1 (x), x e (0, 1) 

O 

On 



The eigenvalue problem associated with (11.11) is given by 



(1.2) u xx — (A 2 + 2Xa — b)u = 0, x G (0, 1), 

^; (1.3) u(0)=u(l) = 0, 

5h I and has received quite a lot of attention in the literature since the pa- 

pers of Chen et al. [CFNS] and Cox and Zuazua [CZj . In the first of 
these papers the authors derived formally an expression for the asymp- 
totic behaviour of the eigenvalues of (11.21) . (II. 3p in the case of a zero 
potential b, which was later proved rigorously in the second of the above 
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papers. Following this, there were several papers on the subject which, 
among other things, extended the results to non-vanishing b |BRj , and 
showed that it is possible to design damping terms which make the 
spectral abscissa as large as desired |CC] . In [F2j the second author of 
the present paper addressed the inverse problem in arbitrary dimen- 
sion giving necessary conditions for a sequence to be the spectrum of 
an operator of this type in the weakly damped case. As far as we 
are aware, these are the only results for the inverse problem associated 
with (11.21) . (jl.3p . Other results for the n— dimensional problem include, 
for instance, the fact that in that case the decay rate is no longer de- 
termined solely by the spectrum [L], a study of some particular cases 
where the role of geometric optics is considered [ALj . the asymptotic 
behaviour of the spectrum [S] and the study of sign-changing damping 
terms |F1] . 

The purpose of the present paper is twofold. On the one hand, 
we show that problem (ll.2p . (11.31) may be addressed in the same way 
as the classical Sturm-Liouville problem in the sense that, although 
this is not a self-adjoinf problem, the methods used for the former 
problem may be applied here with similar results. This idea was already 
present in both |CFNSj and [CZj . Here we take further advantage of 
this fact to obtain the full asymptotic expansion for the eigenvalues 
of (11. 2p . (II. 3p (Theorem [1]). Based on these similarities, we were also 
led to a (regularized) trace formula in the spirit of that for the Sturm- 
Liouville problem (Theorem [4]). 

On the other hand, the idea behind obtaining further terms in the 
asyptotics was to use this information to address the associated inverse 
spectral problem of finding all damping terms that give a certain spec- 
trum. Our main result along these lines is to show that in the case of 
constant damping there is no other smooth damping term yielding the 
same spectrum (Corollary [2]). Namely, we obtain the criterion for the 
damping term to be constant. Note that this is in contrast with the 
inverse (Dirichlet) Sturm-Liouville problem, where for each admissible 
spectrum there will exist a continuum of potentials giving the same 
spectrum |PTj . In particular this result shows that we should expect 
the inverse problem to be much more rigid in the case of the wave 
equation than it is for the Sturm-Liouville problem. This should be 
understood in the sense that, at least in the case of constant damping, 
it will not be possible to perturb the damping term without disturbing 
the spectrum, as is the case for the potential in the Sturm-Liouville 
problem. 

The plan of the paper is as follows. In the next section we set the 
notation and state the main results of the paper. The proof of the 
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asymptotics of the eigenvalues is done in Sections [3] and HI where in 
the first of these we derive the form of the fundamental solutions of 
equation (11.21) . while in the second we apply a shooting method to 
these solutions to obtain the formula for the eigenvalues as zeros of an 
entire function - the idea is the same as that used in [CZ] . Finally, in 
Section [5] we prove the trace formula. 

2. Notation and results 

It is easy to check that if A is an eigenvalue of the problem (II. 2p . 
(11.31) . then A is also an eigenvalue of the same problem. In view of this 
property, we denote the eigenvalues of this problem by A n , n ^ 0, and 
order them as follows 

. . . ^ Im A_2 ^ I m A_i «C Im Ai ^ Im A2 ^ • • • 

while assuming that A_ n = A n . We also suppose that possible zero 
eigenvalues are A±i = A±2 = . . . = X± p = 0. If p = 0, the problem 
(11.21) . (11.31) has no zero eigenvalues. For any function / = f(x) we 
denote (/) := J* f(x) dx. 

Theorem 1. Suppose a G C m+1 [0, 1], b G C m [0, 1], m ^ 1. The 
eigenvalues of U.S\) . U.3\) have the following asymptotic behaviour as 
n — > ±00; 

m— 1 

(2.1) A n = vmi + c 3 n ~ j + 0{n- m ), 

j=0 

were the Cj 's are numbers which can be determined explicitly. In par- 
ticular, 

(a 2 + b) 



(2.2) c = -(a), Cl 



(2.3) c 2 



27ri 

1 



2tt 2 



2 , / v/ 2 , n , a'(l)-a'(0) 



(a(a 2 + b)) - (a)(a 2 + b) 



A straightforward consequence of the fact that the spectrum deter- 
mines the average as well as the L 2 norm of the damping term (as- 
suming b fixed) is that the spectrum corresponding to the constant 
damping determines this damping uniquely. 

Corollary 2. Assume that a G C 3 [0, 1], A n are the eigenvalues of the 
problem U.ty) , U.3\) . the function b G C 2 [0, 1] is fixed, and the formula 
\2. 1\) gives the asymptotics for these eigenvalues. Then the function 
a(x) is constant, if and only if 

c 2 = 2vrici - (b), 
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in which case a(x) = — cq. 

In the same way, the asymptotic expansion allows us to derive other 
spectral invariants in terms of the damping term a. However, these 
do not have such a simple interpretation as in the case of the above 
constant damping result. 

Corollary 3. Suppose 6 = 0, di(x) = a (x) + a>i(x), i — 1,2, where 
a Q (l — x) = a (x), Oj(l — x) — —ai(x), Hi, a G C 4 [0, 1], and for a = Oi 
the problems U.S\) . U.3\) have the same spectra. Then 

(a?) = (oDj ( a oai) = (a al) 

is valid. 

From Theorem [1] we have that the quantity Re(A n — Cq) behaves as 
0(n~ 2 ) as n — > oo. This means that the series 

oo oo 
^ (A n - C ) = 2 ^ Re (^n - C ) 

n=— oo n=l 

converges. In the following theorem we express the sum of this series in 
terms of the function a. This is in fact the formula for the regularized 
trace. 

Theorem 4. Let a G C 3 [0, 1], b G C 2 [0, 1]. Then the identity 

, . a(0) + a(l) 
2^ ( A « - c o) = g ' a ' 

n — — oo 

3. ASYMPTOTICS FOR THE FUNDAMENTAL SYSTEM 

In this section we obtain the asymptotic expansion for the fundamen- 
tal system of the solutions of the equation (11.21) as A — > oo, A G C. This 
is done by means of the standard technique described in, for instance, 
IH Ch. IV, Sec. 4.2, 4.3], [EE Ch. E, Sec. 3]. 

We begin with the formal construction assuming the asymptotics to 
be of the form 

±Xx±J <f> ± (t,X) dt 

(3.1) u±(x, A) = e o j 
where 

m 

(3.2) 0±(x,A) =^0 4 {±) (x)A- 4 + O(A- m - 1 ), m^l. 

i=0 
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In what follows we assume that a G C m+1 [0, 1], b G C m [0, 1]. 

We substitute the series ( 13.11) . (13. 2p into (II. 2p and equate the coef- 
ficients of the same powers of A. It leads us to a recurrent system of 
equations determining which read as follows: 

(3.3) O ±} = a, 

(3.4) 0( ± ) = _I( ±a ' + a 2 + 6)) 

(3.5) # ) = 4( ± ^ ! ' + £^ ±) ^)' i>2 " 

The main aim of this section is to prove that there exist solutions to 
(11. 21) having the asymptotics (13. II) . (13. 21) . In other words, we are going 
to justify these asymptotics rigorously. We will do this for u + , the case 
of u_ following along similar lines. 
Let us write 

m a; 

Ax+E A-/^ +) (t)dt 

t/ m (x, A) = e 1=0 o 

In view of the assumed smoothness for a and b we conclude that U m G 
C 2 [0, 1]. It is also easy to check that 

XJl - X 2 U m - 2XaU m + bU m = \- m e Xx f m (x, A), x G [0, 1], 
(3 ' 6) U m (0) = l, tC(0) = A + f;^(0)A- i ) 

i=0 

where the function f m satisfies the estimate 

\fm(x,\)\^C m 

uniformly for large A and x G [0, 1] 

We consider first the case Re A ^ 0. Differentiating the function u + 
formally we see that 

m 

<(0, A) = A + + (O, A) = A + J2 <ft\W* + 0(X" m ^ 1 ). 

8=0 

Let 

m 

A (A) = A + ^4 +) (0)A"\ 

and A) be the solution to the Cauchy problem for the equation ♦ 

(11.21) subject to the initial conditions 

«+(0,A) = l, u' + (0,\)=A o (X). 
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We introduce one more function w m (x, A) = u + (x, X)/U m (x, A). This 
function solves the Cauchy problem 

{U 2 m w'J + \- m U m e Xx f m w m = 0, x G [0, 1], 

w m (0,A) = l, <(0,A) = 0. 
The last problem is equivalent to the integral equation 
w m (x, A) + \~ m (K m (\)w m )(x, A) = 1, 

X t\ 

(K m (X)w m )(x,X) := / f/" 2 (ti) / U m (t 2 )e xt2 f m (t2,X)w m (t 2 ,X)dt 2 dt 1 . 



Since Re A ^ for ^ t 2 ^ ti ^ 1, the estimate 

\U- 2 (t 1 ,X)U m (t 2 ,X)e xu \^C m 

holds true, where the constant C m is independent of A, t±, t 2 . Hence, 
the integral operator K m : C[0, 1] — > C[0, 1] is bounded uniformly in A 
large enough, Re A ^ 0. Employing this fact, we conclude that 

w m (x) = 1 + C(A" m ), A -> oo, Re A > 0, 

in the C 2 [0, l]-norm. Hence, the formula (13.1 1) , where 

m— 1 

(3.7) + (x, A) = + °( X ~ m ), 

8=0 

gives the asymptotic expansion for the solution of the Cauchy problem 
(fL2l) . (13T6D as A -> oo, Re A ^ 0. 

Suppose now that Re A ^ 0. Let Ai(A), A 2 (A) be functions having 
the asymptotic expansions 

i 

m „ m 

A X (A) = A + X>-* / 4 +) (^)dx, A 2 (A) = A + £ ^ (+) (l)A- 

i=0 q i=0 

We define the function u + (x, A) as the solution to the Cauchy problem 
for equation (11.21) subject to the initial conditions 

« + (l,A) = e Al W u' + (l,X) = A 2 (X)e A ^. 

In a way analogous to the arguments given above, it is possible to 
check that the function u + has the asymptotic expansion (13.11) in the 
C 2 [0,l]-norm as A -> +oo, Re A ^ 0. Hence, w+(0,A) = 1 + 0(A" m ) 
for each m ^ 1. In view of this identity we conclude that the function 
u + (x, A) := m + (x, A)/tt+(0, A) is a solution to (II. 2p . satisfies the condi- 
tion m+(0, A) = 1, and has the asymptotic expansion (13.11) . where the 
asymptotics for + is given in (13 .7p . 
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For convenience we summarize the obtained results in 

Lemma 3.1. Let a G C m+1 [0, 1], b G C m [0, 1]. There exist two linear 
independent solutions to the equation Al.Sfy satisfying the initial con- 
dition u±(0, A) = 1 and having the asymptotic expansions A3. 1\) in the 
C 2 [0, l]-norm as A — > oo ; A G C, where 

m—l 
i=0 

4. ASYMPTOTICS OF THE EIGENVALUES 

This section is devoted to the proof of Theorem [T] and Corollaries |2] 
andEl We assume that a G C m+1 [0, 1], b G C m [0, 1], m ^ 1. 

Let u = u(x, A) be the solution to f 1 1.2 1) subject to the initial condi- 
tions m(0, A) = 0, u'(0, A) = 1. Denote 7o(A) := u(l,A). The function 
7o is entire, and its zeros coincide with the eigenvalues of the problem 
(11.21) . (11.31) . It follows from Lemma 13.11 that, for A large enough the 
function u(x, A) can be expressed in terms of u± by 

_ u+(x,X) -u_(x, A) 
U[X ' >~ <(0,A)-u'_(0,A)' 

The denominator is non-zero, since due to (13. ip 

u' + {0, A) - u'_(0, A) = 2A + 2(a) + C^A" 1 ), A -»• oo. 

Thus, for A large enough 

u+(l,\) A) 
M + (U, A) — tt_(L), A) 

Lemma 4.1. For n large enough, the set 

Q := {A : | Re A| < Tin + 7r/2, | Im A| < nn + 7r/2} 

contains exactly 2n eigenvalues of the problem A3. A3. 3\) . 

Proof. Let 71(A) := 7o(A)e A+ ^. The zeros of 71 are those of 70(A). For 
A large enough we represent the function 71(A) as 

e 2(A+a(0)) _ l 

7i(A) = 72 (A)+ 73(A), 72 := 2(A + a(0)r 

0+(O, A) + 0„(O, A) + 2(1 + A~ 1 a(0))(l - e A_1 ^+(-' A )>; 



73(A) = - 72(A) 



2A(A + a(0)) + <f>+(0, A) + 0_(O, A) 

+ 



e A- 1 (^ f (-,A)) _ e -A- 1 (0_(.,A)) 



2(A + a(0)) + A^(<M0, A) + 0_(O, A)) 
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0±O, A) :=A -1 (0 ± (z, A) - a(x)). 

It is clear that for A large enough the function 73(A) satisfies an uniform 
in A estimate 

| 73 (A)KC|Ar 2 (| 72 (A)| + l). 

One can also check easily that 

| 72 (A)|^C|A|, XedK, 

if n is large enough. These two last estimates imply that (73(A)! ^ 
(72(A) | as A G OK, if n is large enough. By Rouche theorem we conclude 
that for such n the function 71 has the same amount of zeros inside 
Q as the function 72 does. Since the zeros of the latter are given by 
imi — (a), n 7^ 0, this completes the proof. □ 

Proof of Theorem^ Assume first that a G C 2 [0, 1], b G C 1 ! ; 1 ]- As 
was mentioned above, the eigenvalues of problem (11.21) . (II. 3p are the 
zeros of the function 70(A) = 0. It follows from Lemma [4.11 that these 
eigenvalues tend to infinity as n — > 00. By Lemma 13.11 for A large 
enough the equation 70(A) = becomes 

e 2A+<<M-,A)+<M-,A)) _ 
which may be rewritten as 

(4.2) 2A + (</>+(•, A) + </>_(•, A)) — 27mi, n G Z. 

If we now replace <j)± by the leading terms of their asymptotic expan- 
sions we obtain 

(4.3) 2A + 2(a) + ^(A" 1 ) = 27rni, 

A = 7mi — (a) + o(l), n — > 00. 

Hence, the eigenvalues behave as A ~ nni — (a) for large n. Moreover, 
it follows from Lemma 14.11 that it is exactly the eigenvalue A n which 
behaves as 

A n = imi — (a) + o(l), n —> 00. 

It follows from this identity and (14. 3 j) that 

A n = imi — (a) + (9(n _1 ), n — > 00, 

and we complete the proof in the case m — 1. If m — 2, we substitute 
the above identity and (13. ip into (14 .2p and get 

An + (a) + ^(<M +) + 0S" } ) + 0(\- 2 ) = imi, 

X n = imi - (a) - ^L-±$l ] 1 + (n- 2 ). 

imi 
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The last formula and the identities (13.41) yield formulas (12.21) for Co and 
c\. Repeating the described procedure one can easily check that the 
asymptotics (IXTjl . (F2T2|) hold true. □ 

Proof of Corollary^ The coefficients c , C\ in the asymptotics (12.11) 
are determined by the formulas (I2.2p and, by the Cauchy-Schwarz in- 
equality, we thus obtain 

c 2 = (a) 2 ^ (a 2 ) = 2ixi Cl - (b), 

with equality if and only if a(x) is a constant function. This fact 
completes the proof. □ 

Proof of Corollary\E It follows from (E2D, (EJ) that (a 2 ) = (a 2 ), (af) = 
(a^). Now we check that 

(a?) = («o) + (5?>, <a?> = («o) + 3(a a 2 >, z = 1, 2, 
and arrive at the statement of the theorem. □ 

5. Regularized trace formulas 

In this section we prove Theorem HI We follow the idea employed in 
the proof of the similar trace formula for the Sturm-Liouville operators 
in [ESI Ch. I, Sec. 14]. 

We begin by defining the function 

n O-rX 1 -^)- 

n=p+l V A n/ \ A n / 

The above product converges, since 

! (i-A) = i + A2 ~,f,f eAn , 

An/ V A n / |A n | 

and by Theorem [T] we have 

|AJ 2 = vrV - 2vrici + c 2 , + C(n" 2 ), 

(5 1) 

ReA n = c + O(n- 2 ) 

as n — > +oo. Proceeding in the same way as in the formulas (14.8), 
(14.9) in [LSJ Ch. I, Sec. 14], we obtain 

C * (A) sinhA 



$(A) 



A 



OO y- 

n=l ^ 



7r 2 rr - |A n | 2 + 2AReA n 
7r 2 n 2 + A 2 
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Co := (™) 2 *> J] 

ra=p+l 



00 2 2 



I A, 



In what follows we assume that A is real, positive and large. In the 
same way as in \LS\ Ch. I, Sec. 14] it is possible to derive the formula 

(52) ^-tl± ^ : ^ . 

k=l n=l x ' 

Our aim is to study the asymptotic behaviour of ln\l/(A) as A — > +00. 
Employing the same arguments as in the proof of Lemma 14.1 and in 
the equation (14.11) in [LS| Ch. I, Sec. 14], we arrive at the estimate 

k oo ., 



E 

71=1 



ix 2 n 2 - |A TO | + 2A Re A r 
ir 2 n 2 + A 2 

+00 

dt 



C c 



{ (7T 2 n 2 + A 2 ) 



+00 



^ c k X k 



dz 



-k+l 



(ttH 2 + X 2 ) k X k J (tt 2 z 2 + l) fc A 



00 00 / 

^ E^E 

k=3 n=l K 



n 2 n 2 - \X n \ 2 + 2ARe A n 
7r 2 n 2 + A 2 



C(A" 



A — > +00, 



where c is a constant independent of k and n. Let us analyze the 
asymptotic behaviour of the first two terms in the series (15. 2p . As 
k = 1, we have 



00 99 
ir z n z 



E 



|A n | 2 + 2A Re A n , ^ 7T 2 n 2 — \X n \ 2 — 1ix\c\ + Cq 



7r 2 n 2 + A 2 / -— ' 7r 2 n 2 + A 2 

71=1 71 = 1 



+ (27ri Cl - c 2 + 2Ac ) £ 



(5.4) 



71=1 



7r 2 n 2 + A" 



+ 2X^S-2X^y ^ n Y: Xn ~ C0) 
^ ii 2 n 2 + X 2 



71 = 1 



where S := ^^(ReA n — c ). 



n=l 



Taking into account (15.11) . we have 



E 



n=l 



ir 2 n 2 — |A r 



27rici + q 



7r 2 n 2 + A" 



^E 



^ n 2 (ir 2 n 2 + A 2 ) 



7T 3 + A — 3 coth A 



6 



A 



^ CA" 2 , 



EIGENVALUE ASYMPTOTICS FOR THE DAMPED WAVE EQUATION 11 

(Re X n — c )vr 2 n 2 



E 

71=1 



ir 2 n 2 



+ A 2 



71=1 



where the constant C is independent of A. Here we have also used the 
formula 

1 AcothA-1 A^-A" 2 



(5.5) £ 



71=1 



7r 2 n 2 + A 2 



2A 2 



as A — ► +oo. We employ this formula to calculate the remaining terms 

in (15.41) and arrive at the identity 

(5.6) 



7r 2 n 2 — |A n | 2 + 2AReA r 



E 

71 = 1 

as A — > +oo. For k = 2 we proceed in the similar way, 



7r 2 n 2 + A 2 



c +( 2S - c - ^ + iTrci ) A^+C>(A- 2, 



E 

77=1 



vr 2 n 2 - |A„| 2 + 2ARe A n 
7r 2 n 2 + A 2 

00 ' .2^2 



E 



(ttV - |A n | 2 ) 2 
J (tt% 2 + A 2 ) 2 



(7r 2 n 2 - |A ra | 2 )ReA n 2 2 ^ 1 

(tt% 2 + A 2 ) 2 °^(^ 2 + A 2 ) 2 



71=1 
OO 



+ 4A 2 E 



(ReA n ) 2 -c 2 



^ (ir 2 n 2 + A 2 ) 2 ' 

71=1 V ' 

By differentiating (15.51) we obtain 

1 AcothA-2-A 2 (l-coth 2 A) 
^ (vr 2 n 2 + A 2 ) 2 ~ 4A 1 ' 

71 = 1 V ' 

This identity and ( 15. ip yield that as A — > +oo 



E 



(ttV - |A n | 2 ) 2 



^ (ir 2 n 2 + A 2 )' 



^E 



^ (ir 2 n 2 + A 2 V 

77 = 1 V ' 



E 



n=l 

OO 



E 



(vr 2 n 2 + A 2 ) 2 
(ReA n ) 2 - q 



^E 



< (tt% 2 + x 2 y 



^E 



£ (tt% 2 + A 2 ) 
1 



^ CA -3 , 



(7r 2 n 2 + x 2 y 



< c\~ 6 . 



Hence, 



OO / 

(") £ ( 

71=1 V 



7r 2 n 2 - |A w | 2 + 2AReA r 
7r 2 n 2 + A 2 



4\- L + o{x 
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as A -> +00. It follows from (OIL (1531) . (TO . (15771) that 

M(A) = -c - (25 - c + ivrc^A- 1 + 0(A" 2 ), 

(5.8) $( A ) = gog^sinh A _ ^ _ ^ + + o(A _ 2) ^ 

A 

as A — > +00. It follows from (14.11) and Lemma 13.11 that for A large 
enough the estimate 

| 7 o(A)KC|Ar 1 el A l 

holds true. Hence, the order of the entire function 70(A) is one. In view 

00 

of Theorem [1] we also conclude that the series ^2 | A„, | — 2 converges 

n=p+l 

and therefore the genus of the canonical product associated with 70 is 
one. We apply Hadamard's theorem (see, for instance, \Le\ Ch. I, Sec. 
10, Th. 13]) and obtain that 



oc 



7o (A) =e p W$(A), P(X)= ai X + a + 2 ^ |A n |- 2 ReA n , 

n=p+l 

where «i, a are some numbers. Hence, due to (15.81) . it follows that 70 
behaves as 

r< -P(A) \ 

7o(A) = °° e A SmhA [1 - (25 - c + ivrc^A" 1 + 0(A" 2 )] , 
as A — > +00. On the other hand, Lemma [3.11 and (14. II) imply that 

70(A) = ^ [1 + «<M +) > - a(0))A~ 1 + G(A" 2 )] + 0(A^e" A ), 
as A — > +00. Comparing the last two identities yields ot\ = 0, 

00 

and 

-(25 - c + ivrci) = (0i +) > - a(0). 
It now follows from (IO> that 

£ (A n , - c ) = 25 = c + a(0) - (0S +) ) - ivr Cl = Q(0) ± - (a), 



n— — 00 
n#0 



completing the proof of Theorem HI 
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